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Abstract 

A coordinate independent derivation of the Eulerian and Lagrangian strain 
tensors of finite deformation theory is given based on the parallel propagator, 
the world function, and the displacement vector field as a three-point tensor. 
The derivation explicitly shows that the Eulerian and Lagrangian strain ten- 
sors are two-point tensors, each a function of both the spatial and material 
coordinates. The Eulerian strain is a two-point tensor that transforms as a 
second rank tensor under transformation of spatial coordinates and transforms 
as a scalar under transformation of the material coordinates. The Lagrangian 
strain is a two-point tensor that transforms as scalar under transformation of 
spatial coordinates and transforms as a second rank tensor under transforma- 
tion of the material coordinates. These transformation properties are needed 
when transforming the strain tensors from one frame of reference to another 
moving frame. 
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I. BACKGROUND 



The U.S. Army is developing an electromagnetic gun (EMG) for battlefield applications. 
During the past few years, on a recurring basis, Dr. John Lyons (former ARL Director) 
and Dr. W. C. McCorkle (Director of U. S. Army Aviation and Missile Command) have 
requested that I look at some of the physics of the EMG. In the most recent request, I was 
asked to look at stresses in a rotating cylinder. For the case of an elastic cylinder, this is 
a classic problem that is solved in many texts on linear elasticity fl^^BJELH . However, 
when these derivations are examined closely, one finds certain shortcomings J7|. Therefore, 
I spent some time looking at the problem of stresses in elastic rotating cylinders. That work 
resulted in a manuscript @. In the course of this work 0, I had to clearly understand the 
transformation properties of the Largangian and Eulerian strain tensors of finite deformation 
theory. I was quite dissatisfied with the standard derivations of the Largangian and Eulerian 
strain tensors because these derivations take either of two (both unpalatable) approaches. 
In the first approach, shifter tensors are used, which are often defined as inner products 
between two basis vectors at two different spatial locations In this approach, basis 

vectors are not parallel transported to the same spatial location before the inner product 
is carried out. This is unpalatable, even in Euclidean space, unless one is using Cartesian 
coordinates. In the second approach, convected (moving) coordinates are used, and vectors 
and tensors are associated with a given coordinate in the convected (moving) coordinate 
system, rather than being associated with a point in the underlying space. 

In the derivation that I present below, I avoid both of the unpalatable features mentioned 
above. I provide a coordinate independent derivation of the Lagrangian and Eulerian strain 
tensors based on standard concepts in differential geometry: the parallel propagator, the 
world function, and the displacement vector field as a three-point tensor. 

The derivation that I present below is also useful for gaining a basic understanding of the 
role of the unstrained state, or reference configuration, in the definition of the strain tensors. 
Having a firm conceptual grasp of the role of the unstrained state in the definition of the 
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strain tensors is imperative for understanding the behaviour of pre-stressed materials under 
finite deformations in high-stress applications, such as, for example, in the electromagnetic 
rail gun [|16j . 

II. INTRODUCTION 

The theory of stresses in rotating cylinders and disks is of great importance in practical 
applications such as rotating machinery, turbines and generators, and wherever large rota- 
tional speeds are used. In a previous work @, I gave a detailed treatment of stresses in a 
rotating elastic cylinder. This is a classic problem that is treated in many texts on linear 
elasticity theory |2|||,|]|)]||. These treatments linearize the strain tensor in gradient of the 
displacement field, assuming that these (dimensionless) gradients are small. I point out in 
Ref |7j that for large angles of rotation the quadratic terms (in displacement gradient in the 
definition of strain) are as large as the linear terms, and consequently, these quadratic terms 
cannot be dropped. In Ref [0, I provide an alternative derivation of stresses in an elastic 
cylinder that relies on transforming the problem from an inertial frame (where Newton's 
second law is valid) to the co-rotating frame of the cylinder-where the displacements gradi- 
ents are small. During the course of that solution, I had to transform the Lagrangian and 
Eulerian strain tensors of finite elasticity to the (non-inertial) co-rotating frame of reference 
of the cylinder, which is a moving, accelerated frame. This work required the detailed un- 
derstanding of the transformation properties of the Lagrangian and Eulerian strain tensors. 

The standard derivation of these strain tensors is done with the help of shifter ten- 
sors 011]. Shifter tensors are often defined in terms of inner products of basis vectors that 
are located at two different spatial points |8|||. For me, inner products between vectors 
at two different points is an unpalatable operation, even in Euclidean space. In order to 
compute the inner product between two vectors, the vectors must first be parallel trans- 
ported to the same spatial point (unless we are using Cartesian coordinates, in which case 
the derivation becomes coordinate specific). 
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In other treatments, where shifter tensors are not employed in the derivation of strain 



tensors, convected (moving) coordinates are used, see for example fllQ,II|, l2||ni| . When 
using convected (moving) coordinates, the coordinates of the initial undeformed point and 
the deformed point are the same, but the basis vectors change during deformation. In 
derivations of strain tensors using convected coordinates, vectors and tensors are associated 
with a given point in the convected (moving) coordinate system, rather than being associated 
with a point in the underlying (inertial) space. Tensors are absolute geometric objects, and 
they should properly be associated with a point in the underlying space, and not a given 
coordinate, e.g., in moving coordinates. 

In this work, I avoid the unpalatable features of the strain tensor derivation mentioned 
in the above two paragraphs. I derive the strain tensors using the concept of absolute 
tensors, where a tensor is associated with a point in the space-rather than the coordinates 
in a given (moving) coordinate system. I provide a coordinate independent derivation of 
the Lagrangian and Eulerian strain tensors, where I keep track of the positions of the basis 
vectors. The derivation necessarily uses two-point (and three-point) tensors []^|9 yil , |I5 . The 



derivation is based on standard concepts in differential geometry: the parallel propagator 
(a two-point tensor), the world function (a two-point scalar), and the displacement vector 
field (a three-point tensor). This derivation makes clear the transformation properties of 
the strain tensors under coordinate transformations from one frame of reference to a second 
frame that is moving and accelerated (with respect to the first frame). 

The derivation below of the Eulerian and Lagrangian strain tensors makes the trans- 
formation properties (e.g., to a moving frame) clear. Furthermore, this derivation makes 
the role of the reference (unstrained) configuration more clear in the definition of the strain 
tensors. Clarifying this role is of importance for applying finite deformation theory to pre- 
stresssed materials, which are capable of withstanding higher-stress applications, such as in 
rotating machinery |{7,n|. Finally, the derivation presented here allows the generalization 
of the definition of strain tensors to the realm where general relativity applies • 
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III. GEOMETRIC BACKGROUND 



In Euclidean space, a vector can be trivially parallel propagated in the sense that after 
a round trip the vector still points in the same direction. In Riemannian space, the parallel 
displaced vector is not equal to itself after the round trip parallel displacement. In this sense, 
in Euclidean space we need not distinguish the position of a vector because "it always points 
in the same direction under parallel displacement" -even though its components may be dif- 
ferent from point to point because the basis vectors, onto which we project the vector, point 
in a different direction from point to point. So, in Euclidean space the parallel displaced 
(physical) vector (a geometric object) is thought to point in the same physical direction. In 
Riemannian space, however, the situation is quite different. In Riemannian space, when a 
vector is parallel displaced it will in general point in a different direction. The physical test 
is to parallel displace the vector along a curve that returns to the starting point. If there is 
non-zero curvature, as measured by the Riemann curvature tensor, then upon returning to 
its starting point the vector components will be different than the initial vector components 
at the starting point. So, in Reimannian space, it is imperative to specify the position of a 
vector. In Euclidean space, appropriate to material deformations, I also keep track of the 
position of a vector. This additional care in Euclidean space, together with the transforma- 
tion properties of the world function, leads to a clearer understanding of the transformation 
properties of the Lagrangean and Eulerian strain tensors, under transformations from one 
system of coordinates to another that is in relative motion (a moving frame). 

In this section I briefly review the fundamental geometric quantities that naturally arise 
in discussion of deformation, but which are not usually discussed in this context. These 
quantities are the world function (or fundamental two-point scalar of the the space), the 
parallel propagator, and the position vector. This section will also serve to define my nota- 
tion. Each of the quantities mentioned are examples of a class of geometric object objects 
known as two-point tensors, which occur naturally in the discussion of deformations. I have 



found useful discussions of general tensor calculus in Synge and Schild |19] and Synge |15 
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and discussions oriented toward deformation theory in the Appendix by Ericksen in Treus- 
dell and Toupin ||14|| , and in Narasimhan ||, Eringen ||, and Eringen |2T) |. In particular, 



discussion of two-point tensors can be found in Synge ||15|| , Ericksen JTJj], Narasimhan |9| 
and Eringen ||. 

A. The World Function 

The world function was initially introduced into tensor calculus by Ruse pl|,|22|, 



Synge [p3|l , Yano and Muto f24|| , and Schouten p5 |. It was further developed and extensively 



used by Synge in applications to problems dealing with measurement theory in general rela- 
tivity |15| . An application of the world function to problems of navigation and time transfer 
can be found in Ref. [^6|. Compared to the enormous attention given to tensors, the world 



function has been used very little by physicists. Yet, when geometry plays a central role, 
such as in deformation theory, the world function is helpful to clarify and unify the under- 
lying geometric concepts. The world function is simply one-half the square of the distance 
between two points in the space. In applications to relativity and 4-dimensional space-time, 



the space-time is often taken as a general (curved) pseudo-Riemannian space |15|] . In appli- 
cations to deformation of materials, we are concerned with a Euclidean three-dimensional 
space. However, for understanding the transformation properties of displacement vectors and 
strain tensors, it is helpful to use the concept of world function in a Euclidean 3-dimensional 
space described by curvilinear coordinates af , i — 1,2, 3, with a metric gij, which in general 
is a function of position. 

Consider two points in a general Riemannian space, Pi and connected by a unique 
geodesic path (a straight line in Euclidean space) T, given by x % {u), i = 1,2,3, where 
u\ < u < 112, and x l (u) are curvilinear coordinates of the path. The coordinates of point 
Pi = {x\} and point P 2 = {x l 2 }. In general, a geodesic is defined by a class of special 
parameters u', u ■ ■ ■, that are related to one another by linear transformations u' = au + b, 
where a and b are constants. Here, u is a particular parameter from the class of special 



6 



parameters that define the geodesic T, and x l {u) satisfy the geodesic equations 

^+1**^ = (1) 
du 2 3 du du 

Using Cartesian coordinates z k (rather than general curvilinear coordinates x k ) in Euclidean 
space, the Christoffel symbol Y^ k = 0, and the solution of Eq. ([[]) is simply 

z\u) = z? + (z« - *?) (2) 

u 2 - U\ 

where u\ < u < u 2 , % = 1, 2, 3 and the Cartesian coordinates of points Pi and P2 are zf and 
z% , respectively. In a general Riemannin space, the world function between point Pi and P2 
is defined as the integral along V in arbitrary curvilinear coordinates x % by 

n(p 1 ,ft) = I(« 2 -« 1 )£ a ^^d« (3) 

The value of the world function has a simple geometric meaning: it is one-half the distance 
between points Pi and P%. Its value depends only on the eight coordinates of the points 
Pi and P2. The value of the world function in Eq. (|3|) is independent of the particular 
special parameter u in the sense that under a transformation from one special parameter u 
to another, u', given by u = au' + b, with x l (u) = x l (u(u')), the world function definition in 
Eq. (|3p has the same form (with u replaced by u'). 

The world function is the fundamental two-point invariant that characterizes the space. 
It is invariant under independent transformation of coordinates at Pi and at P2. For a given 
space, the world function between points Pi and P2 has the same value independent of the 
coordinates used, which makes it a useful coordinate independent quantity. In Euclidean 
space, using Cartesian coordinates, the world function has the simple form 

Q(zi4) = ^S ij Az i Az j (4) 

where dij is the Euclidean metric with only non-zero diagonal components (+1, +1, +1), 
and Az l = {z\ — z\), i = 1,2,3, where z\ and z % 2 are the Cartesian coordinates of points 
Pi and P2, respectively. (I use the convention that all repeated indices are summed, unless 
otherwise stated.) 



The world function has a number of interesting properties, see Synge fll5| . Calcula- 
tions of the world function for spaces other than Euclidean spaces, namely four- dimensional 
space-time, can be found in Refs. |15| , |26|p7| , |28|p9| . In what follows, I restrict myself to a 
three-dimensional space. By transforming to a new system of coordinates, say spherical 
coordinates, 



x = r cos 9 cos< 
y — r cos 6 sin ( 
z = r cos 6 



(5) 
(6) 
(7) 



the world function in Eq. (|]) can be expressed as a function of spherical coordinates of point 
Pi = (r 1 ,9 1 ,4> 1 ), and P 2 = (r 2 , 02, fa)- 

Consider a geodesic given by Eq. ([]]) in a general 3-dimensional Riemannian space. The 
covariant derivatives of the world function have two important properties: 

dn(Px,p 2 ) 



dx 2 



(u2 - ui) {^ij ^-^J = L X h 

-(u 2 - Ul ) {m d ±) p j 



(8) 
(9) 



where 



L=[2Q(P 1 ,P 2 )] 1/2 = r ds= ( 

J Pi Ju 



I Ul 



9ij 



dx l dxi 
du du 



nl/2 



du 



(10) 



is the length of the geodesic between P\ and P 2 , gij is the metric in coordinates x l , and 
Aj x and \ 2 are components of the unit tangent vectors at end points Pi and P 2 (assuming 
non-null geodesies |fL5|j): 



/ dx^ 



A; 



9ij 



dx j \ 



(11) 
(12) 



\ ds I D 

where the relation between parameter u and arc length s is given by Eq. (|T0"D. In Eq. (^|) 
and (^|), the covariant partial derivatives with respect to x\ and x l 2 are done with respect to 
the coordinates of points Pi and P 2 . See Fig. ([[]). 
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For the special case of interest in deformation of materials, the space is three-dimensional 
Euclidean, the geodesic is a straight line, and the vectors and X i2 are colinear, although 
I still consider them as existing at distinct points. Using Cartesian coordinates and the 
explicit form of the world function in Eq. (f|), Eq. (|8|) and (|9[) take the form 

dQ(P u P 2 ) 



d z. 



2 



Vt i2 = 5 i:j (j4 - 4) = LX i2 (13) 



an( Pl ,p 2 ) = ^ iUi = _ 5i . u _ z{ \ = _ LXii (14) 



dz\ 

where I used Synge's short-hand notation for the components of the covariant partial deriva- 
tives by putting subscripts on the indices to indicate which coordinates were differentiated. 
This short-hand notation is particularly convenient to show the transformation properties of 
the world function and to indicate the spatial location of vectors and tensors. For example, 
the quantity Qi 2 transforms as a vector under coordinate transformations at P 2 and as a 
scalar under coordinate transformations at point Pi. The quantity \ 2 is a vector located at 
point P 2 . Note that by virtue of their definitions in the left side of Eq. (|H|) and (|I4]), the 
right sides are two-point tensors, whose components are functions of coordinates at point 



Pi and P 2 . For example, the right side of Eq. (13) is a product of a two-point scalar L, and 
a one point vector, A^ at Pi. 



B. Parallel Propagator 

Given a vector with components v %1 at point Pi, the vector is said to be parallel prop- 
agated from Pi to P 2 along a geodesic curve C specified by x l (u), U\ < u < u 2 , where 
Pi = {x l (ui)} and P 2 = {x*(u 2 )}, when its covariant derivative is zero along this curve: 

dv l • ■ dx k 

— + r s y — = o 15 

du Jk du v ' 

Equation (|15|) is a mapping: given the components of a vector, v H at point Pi, we obtain 
the components v 12 of the parallel transported vector at point P 2 by solving Eq. flT5|). It 
is convenient to define a two-point tensor, g 12 ^, called the parallel propagator fl5|], which 



gives the components of a vector under parallel translation of the vector from point P\ to 
point P 2 . Given a vector with components v 11 at point Pi, the propagator g 12 ^ relates the 
components of this vector at Pi to the components v 12 of this same vector after parallel 
translation to point P 2 

v»=g\v il (16) 

In a general Riemannian space, the components of the vector at point P 2 depend on the 
path of parallel translation from Pi to P 2 , in the sense that the path must be a geodesic by 
the definition of the parallel propagator. However, in Euclidean space these components are 
completely path independent; the components depend only on the end points Pi and P 2 . 

A vector is considered as a geometric object, which means that it is independent of 
coordinate system. In a Riemannian space, under the operation of parallel propagation a 
vector changes in such a way that its magnitude stays the same but its absolute direction can 
change because of the curvature of the space |!0| . The direction of the parallel propagated 
vector is of course referred to the local basis vectors. That the vector direction changes under 
parallel translation can be understood by taking a vector at point P and parallel translating 
it over a curve that returns to point P. When compared at point P, the components of the 
initial vector and the round-trip-parallel-transported vector will (in general) be different. It 
is in this sense that a vector changes its direction under parallel transport. 

As mentioned above, the change in the vector that results under parallel transport de- 
pends on the path of parallel propagation (a geodesic). Two vectors that are parallel prop- 
agated along the same path will maintain the angle between them along the path. 

In a Euclidean space, a vector (the geometric object) is considered to be unchanged when 
parallel propagated. The only thing that happens is that the components of the vector on 
the local basis must change according to what is required to keep the vector "pointing in 
the in same direction" . 

In Euclidean space, the parallel propagator in Cartesian coordinates is trivial-its com- 
ponents are just the components of a delta function. The components of a vector at point 
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Pi are related to the components of the same vector parallel translated to point P2 by the 
propagator (whose components are given in a Cartesian coordinate basis): 

+1 i=j 

*\ = ( 17 ) 

Equation ([17]) agrees with our notion from elementary geometry that in Cartesian coordi- 
nates the vector components are constant under parallel propagation. However, using the 
parallel propagator in Cartesian coordinates, I can, for example, compute the propagator 
g 12 j 1 in curvilinear coordinates x l = (r,9,(p) given in Eq. by the two-point tensor 

transformation rule 



dx^dz^P,) 

9 n--d^r—^r 5 - ^ 



The parallel propagator g %2 ^ is a two-point tensor because it transforms as a vector under 
coordinate transformation at point Pi and under coordinate transformation at point Pj. 

In Cartesian coordinates, when the points are made to coincide, P 2 — > Pi, the propagator 
reduces to a Kronecker delta at point Pi: 8™% — > 5 m n (Pi). In general curvilinear coordi- 
nates, when the points Pi and P2 coincide, the parallel propagator reduces to the mixed 
components of the mertic tensor g 12 ^: 

hia g i2 h -> ^(Pi) (metric at P x ) (19) 

The mixed components of the metric tensor at P 1; g l ■ (Pi) = g lk (Pi) g^j (Pi), are a Kronecker 
delta-a unit tensor whose components are the same in all systems of coordinates. Indices 
can be lowered on two-point tensors using the appropriate metric. For example, the index i 
of the propagator g l2 j i can be lowered by using the metric tensor at point P 2 : 

When the points are made to coincide, P2 — ► Pi, the covariant components of the propagator 
becomes the covariant components of the metric tensor at Pi, g^n ~* 9kj(Pi), where gkj{Pi) 
is the metric at Pi. 
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The covariant derivatives of the world function Q(Pi,P2) between points P\ and P2 are 
related to the parallel propagator by |T5f| 



^hh = 9hh (metric at Pi) (21) 
Q hh = tt j2il = -g ilj2 = -g j2il (parallel propagator) (22) 
fli 2 j 2 = Qi 2 j 2 (metric at P 2 ) (23) 



Other useful properties of the parallel propagator are discussed by Synge 15 



C. Position Vector 

The position vector occupies a central role in deformation theory. For this reason, I 
discuss it in detail below. In elementary geometry, a point P can be identified by its position 
vector r, which can be specified in Euclidean-space Cartesian coordinates as 



r 



z n i n (24) 



where z n are the Cartesian components of the vector r and also the Cartesian coordinates 
of the point P. In terms of general coordinate basis vectors e n = d/dx n associated with the 
curvilinear coordinates x n , the vector r is given by 

r = z n A m n (P) e m (P) = ( m e m (P) (25) 

The position vector is a geometric object at point P. Among all basis vectors, the Cartesian 
basis vectors i n are unique in that they are usually not associated with a particular spatial 
point. However, when we express these Cartesian basis vectors in terms of curvilinear basis 
vectors e n , then we must imagine that these basis vectors exist at a particular point P. 
Hence, the transformation between the Cartesian basis vector i m and curvilinear coordinate 
basis vectors e m at point P associated with coordinates x % is given by 

i n (P) = A m „(P)e m (P) (26) 

where the matrix A™(P) depends on the coordinates of point P: 
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^n(P) = J^{P) (27) 

In Cartesian coordinates, the components of the vector r are simply the Cartesian co- 
ordinates z n of point P. The three numbers (z l ,z 2 ,z 3 ) transform as the components of a 
vector under orthogonal coordinate transformations. Note that in curvilinear coordinates, 
the components of the position vector, ( m , are not the curvilinear coordinates of point P. 
Also, note that the position vector r of point P has a magnitude equal to the Euclidean 
length from the origin of coordinates, say point O, to point P. The position vector of point 
P is a geometric object at point P, however, it also depends on the point O. This depen- 
dence on point O is coordinate independent. Therefore, the position vector of point P is a 
two-point tensor; it depends on point P and on point O. The transformation properties of 
the position vector are that of a scalar when a change of coordinates is made at point O and 
the transformation is that of a vector when coordinates at point P are changed. 

In a Riemannian (a generalization of Euclidean space) space, the components of the 
position vector r«(P) at point P can be defined in terms of the covariant derivative of the 
world function 

r iP = dQ i°: P) = V ip (0,P) = [2n ip (0,P)] 1/2 r\P) (28) 

Zp 

1 /2 

where r %p is a unit vector at point P and [2 Qi p (0, P)\ is the length of the geodesic from 

1 /2 

point O to point P. For the case of Euclidean space, [2 fl, ip {0, P)\ is the length of the 
straight line OP. Equation ( p8[ ) shows explicitly that the position vector, r ip , is a two-point 
tensor. 

D. Displacement Vector 

Consider an elastic body that undergoes a finite deformation in time. The deformation 
can be specified by a flow function or displacement mapping function 

z k = z k (Z m ,t) (29) 
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where the coordinates z k (here taken to be Cartesian) of a particle at point Q at time 
t are given in terms of the particle's coordinates Z k of point P in some reference state 
(configuration) at time t = t a , so that z k (Z m ,t Q ) = Z k . I assume the deformation mapping 
function has an inverse, which I quote here for later reference 

Z m = Z rn (z k ,t) (30) 

I assume that both the coordinates z k and Z m refer to the same Cartesian coordinate 



system |33 



In deformation theory, the initial position of the particle in the medium at point P is 
specified by a vector 

R(P) = Z m ki k (P) (31) 

and the final position is specified by a position vector 

r(Q)=z k i k (Q) (32) 

where the quantities ifc(-P) and ik{Q) are the Cartesian basis vectors at point P and point 
Q, respectively. Conventionally, the deformation of a medium is described by specifying the 
displacement "vector field", which is defined as a difference of these two position vectors. 
However, the basis vectors h(P) an d h(Q) are at different points in the space. Since vectors 
can be subtracted only if they are at the same point, I must parallel translate i&(P) to point 
Q, or, parallel translate ik{Q) to point P. Depending on which mapping I choose, I arrive 
at the Eulerian or the Lagrangian displacement vector. 

First, I parallel translate vector R(P) to point Q, and then subtract the vectors at point 
Q, see Fig. (|2|). This procedure defines the components of the Eulerian displacement vector 
at point Q: 

u(Q) = r(Q) - R(Q) (33) 

This Eulerian displacement vector in Eq. ( |3~5| ) is often called the displacement vector in 
the spatial representation P,P,|3l|. Alternatively, I can parallel translate the vector r(Q) to 
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point P, and then subtract the vectors at point P. This procedure defines the Lagrangian 
displacement vector at point P: 

U(P) = r(P) - R(P) (34) 



This Lagrangian displacement vector in Eq. (|34| ) is often called the displacement vector in 
the material representation [^|,|9|,|3ll . Equations ( |3~3"|) and (^) show that these two vectors 
are actually referred to basis vectors at different points. In fact, the two vectors u(Q) 
and U(P) are related by parallel translation. In a Euclidean space, these vectors are the 
same geometric objects but they are expressed in terms of basis vectors located at different 
positions 

In order to further clarify the transformation properties of these two displacement vec- 
tors, I use the position vector as discussed in the previous section. Consider the deformation 
mapping function in curvilinear coordinates, x k (X m ,t). This function specifies the coordi- 
nates x k (point Q) of a particle at current time t in terms of the coordinates X k (point P) 
of the particle in the reference configuration at time t = t , so that 

x k {X m ,t )=X k (35) 

In addition, there exists a straight line (a geodesic) T connecting the points P and Q. 

The covariant components of the position vector of point P, R(P) = R n e n (P), in 
curvilinear coordinates x l are given by (see Eq. (p8|)) 



Ri P = 9 ^Q°Jp = n ip (0,P) = [2 n(0,P)] 1/2 R lp (36) 



where Ri p are the components of the unit vector at point P tangent to V that connects 
point P and Q. Similarly, the covariant components of vector r(Q) = r n e n (Q) in curvilinear 
coordinates x % are given by 



n Q = 'J = n iQ (o,Q) = [2n(0,Q)] l/2 f iQ (37) 



lQ ~ dx^ 

where fj Q are the components of the unit vector at point Q tangent to T. From Eq. (|36] 
and ( p7|) it is clear that both P, p and r, Q are two-point tensor objects. The quantity Ri 1 
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depends on point O and P and transforms as a vector under coordinate transformations at 
P and as a scalar under coordinate transformations at point O. The quantity r iq transforms 
as a vector under coordinate transformations at point Q and as a scalar under coordinate 
transformations at point O. 

The components of the Eulerian displacement vector in Eq. ( |33"D (at point Q) are defined 
in terms of the components of R lp parallel translated to point Q: 

BjQ = g ^ jp Bjr (38) 
= g\ Q jp (0,P) (39) 

where fi jp (0, P) = g^ k {P) VLk P {0, P), and gi k (P) is the metric at point P in coordinates x % . 
The contravariant components of the Eulerian displacement vector in Eq. (1331) are given by 



= (0,Q)- g% n» (O, P) (40) 

Similarly, the components of the Lagrangian displacement vector are given by 

U*r = g\ (0,Q)- (O, P) (41) 

The vectors whose components are U lp and u lQ , are related by parallel transport along the 
geodesic V connecting P and Q (not along the particle displacement line given by Eq . (|29"D ) . 
Transporting U lp to point Q 

U lQ = 9% P U kp (42) 

= g\ [g 3P kQ ^ kQ (o, Q) - (o, p)} (43) 

= 8\{Q)n k i{o, Q) - ^ p tf*{0, P) (44) 

= (45) 



where in the transition from Eq. (|43| ) to ( ^4|) I have used the identity satisfied by the parallel 
propagator: 

8\{Q) = g\g j \ Q (46) 
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where S\(Q) is a unit tensor (delta function) at point Q. Equation fl46|) is the statement 
that parallel propagation of a vector has an inverse, so the result that U l (P) and u l (Q) are 
related by parallel transport is true in both Euclidean and Reimannian spaces. 



IV. STRAIN TENSORS 

At time t = t D , consider two particles in the medium that are at positions Pi and P 2 , 

1 /2 

respectively, and are separated by a finite distance AS* = [2f2(Pi, P 2 )] . At a later time 
t > t Q these particles have moved to new positions Qi and Q 2 and are separated by a distance 
As = [2fi(Q 1; Q 2 )] 1/2 , see Fig. ©. 

As a measure of strain, I take the 4-point scalar 

*(P l5 P 2 ; Q x , Q 2 ) = (As) 2 - (AS) 2 = 2Q(Q U Q 2 ) - 2Q(P U P 2 ) (47) 

Note that \l/(Pi, P 2 ; Qi> Q2) depends on four points Pi, P 2 , Qi, and Q 2 , and by virtue of 
its definition in terms of the world function, \1/ is a true four-point scalar invariant under 
separate coordinate transformations at each of these four points. In Cartesian coordinates, 
Eq. ( pETf ) becomes 



9{4, 4; Z{, Z{) = fiy (z\ - z\) (z{ - z{)) - 5 tJ (Z\ - Z\) (Z{ - Z{)) (48) 

where z\ = z 1 (Z] a ,t) and z\ = z t (Z' 2 TL ,t) and they are related to the reference configuration 
at t = t by 

z l (Z™,t ) = Z{ (49) 

and 

z l (Z™,t ) = Zi (50) 

and z l (Z m ,t) is the deformation mapping function in Cartesian coordinates, given in 
Eq. (]2"9"D. If I consider the particle positions Pi and P 2 as separated by an infinitesimal 
distance, then, by assuming continuity in the medium and a finite time t — t Q , the points 
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Qi and Q2 are also infinitesimally separated. However, because t — t is finite, the distance 
between Pi and Qi is finite (not infinitesimal). Expanding Eq. ( |50| ) about the initial position 
of the first particle 

4 = z\Ztt) + §§j(Z?,t) (Zi -Z{)+... (51) 

using z\ = z l (Z±,t), leads to the relation between spatial (Eulerian) coordinates z % and 
material (Lagrangian) coordinates Z k 

Az i = ^ J (Z?,t)AZi + ... (52) 



where Az % = z\ — z\ and AZ % = Z\ — Z\. Using Eq. Qo"2l ) in Eq. (^) leads to the measure 
of strain in Cartesian coordinates 

( M 2 - (AS) 2 = (<% ^ ||r - s™) Azm Azn + ■ ■ ■ ( 53 ) 

= 2 E mn AZ m AZ n + ■■■ (54) 

where the quantity in parenthesis is the Lagrangean strain tensor, E mn . The higher order 
terms in AZ m are small because I can always choose the two initial points P\ and P2 



arbitrarily close together. From Eq. (|53|)-(5_4), it is clear that the Lagrangean strain tensor 



is a two-point tensor, depending on initial point P (in the reference configuration at t = t Q ) 
and point Q (at time t). Note that in Eq. (|54j ) there is no restriction to short times t — t a , 
since I can always choose AZ n sufficiently small. 

The Eullerian tensor can be obtained from Eq. ([3]) by using the fact that the flow 
function in Eq. (|2TS| ) has an inverse. Using 

AZ n = -—-(z\t) Az i (55) 

dz* 

I get the measure of strain in terms of the Eulerian strain tensor e mn : 

(A S ) 2 - (AS) 2 = [5 mn - S kl |^ HI) Az m Az n H (56) 
= 2 e mn Az m Az n + --- (57) 
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A. Strain Tensor Derivation in Curvilinear Coordinates 



I return to the definition of the measure of strain given in Eq. (fI7|). In the reference 
configuration at t — t Q , consider two particles at points Pi and P 2 with curvilinear coordi- 
nates X\ and X 2 . At a later time t, these two particles are at positions Qi and Q2 with 
curvilinear coordinates X\ and £ 2 , respectively. Consider the first term on the right side of 



Eq. (f£TD, which is an integral along a straight line QiQ 2 : 

CXQi, Q2) = 7^2 - w 1 ) / 9ij W W dw (58) 



where U l = dx l (w)/dw and where the geodesic (straight line) is parametrized by x l (w), with 
w± < w < W2 and the end points are given by x\ = x % (w\) and X2 = x l (w2), see Fig. |||. 
The function x l {w) is a solution of the geodesic Eq. ([!]). In the case of Eulidean space, and 
assuming points Pi and P 2 are arbitrarily close, the geodesic in Eq. (|5|) is a straight line 
given by 

xHw) = x\ + w ~ Wl (4 _ x [) (59) 

W 2 - Wi 



The flow function in Eq. ( p9|) maps the points Pi and P2 into the points Qi and Q 2 . The 
points Qi and Q 2 depend on time t. With reasonable continuity assumptions, and the 
straight line given in Eq. (|5TS|) with U lQ = k(x 2 — x\) = k Ax 1 and k = (u>2 — the 
world function in Eq. ( pq) can be approximated by 

M(Qi, Q 2 ) = t:{w2- Wi) 9ij(Qi) kAx 1 kAx 3 / dw (60) 
= ±g ij (Q 1 )Ax i Axi (61) 

Similarly, the second term on the right side of Eq. (|47| ) can be approximated as 

Q(P U P 2 ) = i ^(Pi) AX i AX 3 (62) 

where the coordinates X 1 are the undeformed ones and AX 1 = X % 2 — X\. The measure of 
strain in Eq. (f|7]) is then 
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(As) 2 - (AS) 2 = gij (Q) Ax 1 Ax j - gij (P) AX 1 AX j (63) 



or using the flow function in Eq. (f29|). 

(As) 2 - (AS) 2 = (gaiQ) |^(P, Q) |J(P, Q) - g M {P)^ AX k AX 1 (64) 

Note that x l and X % refer to to the same system of coordinates. I have dropped the subscripts 
on Q and P since Qi and Q2, and Pi and P2, are infinitesimally close, respectively. The 
quantity in parenthesis on the right side of Eq. (EM is twice the Lagrangian strain tensor: 



3t % ot^ 

2E kplp =g ij (Q) MK (P,Q) MJ (P,Q) - g kl (P) (65) 

From Eq. (|6~5|) , it is clear that the Lagrangian strain tensor is a two-point tensor. Under 
transformation of coordinates at point P, Ek P i P is a second rank tensor, while under trans- 
formation of coordinates at point Q, it is a scalar. The deformation gradient, dx' l /dX h , is 
itself a two-point tensor, as can be seen by its transformation property when coordinates at 
P and Q are changed. 

It is interesting to note that the Lagrangian strain tensor Em is conventionally thought to 
be a function of material coordinates, X\ which coincide with the point P (in the reference 
configuration) @,[|. The tensor E kp i P can be taken to be a function of only the material 
coordinates by using the flow mapping function in Eq. fl29|), which provides a mapping 
between all points P and their images, points Q, under the deformation. I do not pursue 
this interpretation below. 

The first term in Eq. (|65|) is the Green deformation tensor: 

Cu = gn{Q)^s(P,Q)^i(P,Q) (66) 

The point P is in the reference configuration at time t = t Q and the point Q is in the 
deformed state at time t. In Eq. (|66"D, the Green deformation tensor is naturally a second 
rank tensor with respect to material coordinate transformations (point P) and it is a scalar 
spatial coordinate transformations (point Q). However, the Green tensor is conventionally 
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taken |||9] as a function of material coordinates by using the flow mapping function in 
Eq. ©. 

Returning to Eq. fl64|), and using the flow mapping function in Eq. (p9|), we can obtain 

9ij(Q) - 9ki(P) ^~(P, Q) ^~(P, Q)J Ax* Ax> (67) 

where the Eulerian strain tensor Ci Q j Q is given by 

dX k dX l 

2e iQiQ = </«(Q) - 9h(P) -g-(P) g-(P) (68) 

The Eulerian strain tensor &i Q j Q is a two-point tensor that is second rank with respect to 
spatial coordinates at point Q and is a scalar with respect to material coordinates at point P 
. Once again, however, using the flow mapping function in Eq. (|29|), e iQ j Q can be imagined 
to depend on the material coordinates x l (point Q) of the deformed state. 

V. STRAIN TENSORS IN TERMS OF THE DISPLACEMENT FIELD 

The Eulerian and Lagrangian strain tensors can be expressed in terms of the displacement 
fields u lQ and U lp in Eq. (^) and (|4ll) . In terms of covariant components, the displacement 
field vector at Q is given by 

u iQ = n iQ (0,Q)-g iQjp Wr(0,P) (69) 

The covariant derivative of Ui Q with respect to point Q (coordinates x h ) is 

dX l 

Ui Q -,k Q = tt iQkQ (O, Q) - 9iQ P n jplp (O, P) (70) 

where I used the chain rule for covariant differentiation since Q(0,P) = Q(0, P(X)) and 
the material coordinates X 1 = X l (x k ) are functions of the spatial coordinates x k . It is 
clear that the chain rule must be used in Eq. ([T(]) by considering Cartesian coordinates. In 
Eq. (]7(]|), I also used the fact that in Euclidean space the parallel propagator is a constant 
under covariant differentiation. The notation Qi Q k Q (0,Q) means the (i, k) component of 
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the second covariant derivative of the world function at point Q. In Euclidean space, these 
second covariant derivatives are simply related to the parallel propagator, see Eq. (p3|) . 
Using Eq. (|23|), the covariant derivative of the displacement field in Eq. (|70| ) becomes 

OX 1 

u iQ -k Q = 9ik{Q) ~ 9 iQ 3P gji{P) -Ty; (71) 

dX l 

= 9ik(Q) - g iQ i P -7^7 (72) 

where in the last line the metric at P, gji(P), was used to lower the index on the propagator. 
Now I multiply Eq. flT2"|) by the parallel propagator g mpl Q ; sum on iQ, and solve for the inverse 
displacement gradient 

F)X m 

= 9 m \ ~ g m *> u lQ , kQ (73) 

From Eq. (|73|), it is clear that in Euclidean space, the deformation gradient dX m /dx k is 
simply related to the covariant derivative of the displacement field, Ui Q -k Q - Note however, 
that in a Riemannian space, for finite deformations, it is generally not possible to solve for 



the deformation gradient [18|. Now, using Eq. ([73|) in Eq. (0), I find the expression relating 
the two-point Eulerian strain tensor &i Q j Q = ^i Q j Q {P,Q) to the covariant derivatives of the 
three-point displacement field 

(As) 2 - (AS) 2 = [u iQ . jQ + u jQ . iQ - g k \Q) u kq . iQ u lQ;jQ ] Ax 1 Ax j (74) 
= 2 e iQjQ Ax 1 Ax j (75) 

where 



ei Q]Q 2 



Ui Q\3Q ^ U 3Q\ i Q 9 (Q) U k Q ;iQ U Iq Jq (76) 



Equation (|76| ) explicitly shows that the Eulerian strain tensor Ci Q j Q is a function of two 
points, material coordinates at point P and spatial coordinates at point Q. Note that Ci Q j Q 
is not a function of point O, since by Eq. (|7T|) the covariant derivative u iQ] j Q does not depend 
on point O. From Eq. (^) it is also clear that Gi Q j Q transforms as a second rank tensor 
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under spatial coordinate transformations and that it transforms as a scalar under material 
coordinate transformations. 

An analogous relation can be obtained for the Lagrangian strain tensor by considering 
the covariant derivative of the displacement field 

Sx^ 

u ip ; kp =gi/ Q n jQlQ (o,Q)-^-n ipkp (o,p) (77) 

Using the relations between the second covariant derivatives of the world function and 
propagator in Eq. (|2l|)-(|23D, and solving for the displacement gradient I get 

Inserting the displacement gradient in Eq. (^) into Eq. ( |B~S[ ) gives an expression for the 
Lagrangian strain tensor E mpnp = E mpnp (P,Q), 

(As) 2 - (AS) 2 = [U mp , np + U np , mp + gV(P) U lp , mp U JP , np ] AX m AX n (79) 
= 2E mpnp AX l AX> (80) 

where 



E — - 

J ->mprip 2 



Urnp;np U np - mp -\- Cj ^ (P) Ui p - mp Uj p ^ np (^1) 



Equation ( [79] ) shows that the Lagrangian strain tensor E mpnp is a function of the material 
coordinates at point P and spatial coordinates at point Q. The Lagrangian strain tensor 
transforms as a scalar under spatial coordinate transformations at point Q and as a second 
rank tensor with respect to material coordinate transformations at point P. Note that there 



are minus sign differences in Eq. ([76]) and (qT|) - Finally, comparing Eq. (|74|) and (f79[), we 
have the well-known relation between the two strain tesnors 

dx l dx j . . 

E mp np = Gi Q j Q (82 J 

Equation ( [52]) provides a complicated relation between the two two-point strain tensors. 
While the displacement fields u %Q and U lp are related to each other by parallel transport, 
see Eq. (pEU|)-([n|), the strain tensors E mpnp and ej QjQ are related by two-point deformation 
gradient tensors, dx l /dX m , in Eq. ([82]). 
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VI. SUMMARY 



Conventionally, the Eulerian and Lagrangian strain tensors are derived either by using 
shifter tensors or by using convected (moving) coordinates. The definition of the shifter 
tensor makes use of a scalar product between vectors at two different points in space (without 
first parallel translating one vector to the position of the other). When convected coodinates 
are used, vectors and tensors are associated with given coordinates in the convected system 
of coordinates, rather than being associated with a gven point in the underlying space. As 
discussed in the introduction, both of these features are undesireable, when we need to 
understand the transformation properties of the strain tensors from an inertial frame to a 
moving frame. These transformation properties are also needed in order to generalize the 
strain tensor to Riemannian geometry for applications to general relativity. 

I have provided a derivation of the Eulerian and Lagrangian strain tensors for finite de- 
formations using the concepts of parallel propagator, the world function of J. L. Synge and 
the three-point displacement vector field. This derivation avoids the undesireable features 
mentioned above. The derivation shows that the Eulerian strain tensor is a two-point object 
that transforms as a scalar under transformation of material coordinates and as a second 
rank tensor under transformation of spatial coordinates. The derivation also shows that the 
Lagrangian strain tensor behaves as a scalar under transformation of spatial coordinates 
and as a second rank tensor under transformation of material coordinates. These trans- 
formation properties are useful in understanding how these strain tensors transform from 
one frame of reference to another moving, non-inertial frame. The formulation presented 
here of the transformation properties of these tensors is also useful for understand the role 
of the reference (unstrained) configurartion in pre-stressed materials, as discussed in the 
introduction. 
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Figure 1 

FIG. 1. A geodesic path is shown in 3-dimensions, with tangent unit vector at the ends. 
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Figure 2 
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FIG. 2. The initial and final position vectors, R(P) and r(Q), respectively, are shown as well 
as their respective parallel translated vectors, R(-P) and r(Q). Also, shown by a solid curve is the 
actual displacement path of a representative particle of the medium, labeled by z k = z k (Z m ,t). 
The dashed straight line is the line (geodesic) connecting the initial and final particle positions. 
The Eulerian displacement vector is u = r(Q) — R(Q), which is the difference of two position 
vectors at point Q. The Lagrangean displacement vector, U = r(P) — R(P), is the difference of 
two position vectors at point P. 
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Figure 3 

FIG. 3. The position Pi and P2 of two particles is shown in the reference configuration at 
t = t , and the positions Q\ and Q2 of the same two particles is shown at later time t > t Q . The 
path of the particles is shown in solid lines and their displacement is shown in dashed lines. 
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